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ABSTRACT 



Aims. Recent studies of the optical/UV and X-ray ephemerides of XI 822-371 have found some discrepancies in the value of the 
orbital period derivative. Because of the importance of this value in constraining the system evolution, we comprehensively analyse 
all the available optical/UV/X eclipse times of this source to investigate the origin of these discrepancies. 

Methods. We collected all previously published X-ray eclipse times from 1977 to 2008, to which we added the eclipse time observed 
by Suzaku in 2006. This point is very important to cover the time gap between the last RXTE eclipse time (taken in 2003) and the 
most recent Chandra eclipse time (taken in 2008). Similarly we collected the optical/UV eclipse arrival times covering the period 
from 1979 to 2006, adding a further eclipse time taken on 1978 and updating previous optical/UV ephemeris. We compared the X-ray 
and the optical/UV ephemeris, and finally derived a new ephemeris of the source by combining the eclipse arrival times in the X-ray 
and optical/UV bands. 

Results. The X-ray eclipse time delays calculated with respect to a constant orbital period model display a clear parabolic trend, 
confirming that the orbital period of this source constantly increases at a rate of P ^ = 1.51(7) x 10~ 10 s/s. Combining the X-ray 
and the optical/UV data sets, we find that P orb = 1.59(9) x 10~ 10 s/s, which is compatible with the X-ray orbital solution. We also 
investigate the possible presence of a delay of the optical/UV eclipse with respect to the X-ray eclipse, finding that this delay may not 
be constant in time. In particular, this variation is compatible with a sinusoidal modulation of the optical/UV eclipse arrival times with 
respect to the long-term parabolic trend. In this case, the optical/UV eclipse should lag the X-ray eclipse and the time-lag oscillate 
| about an average value. 

Conclusions. We confirm that the orbital period derivative is three orders of magnitude larger than expected from conservative mass 
transfer driven by magnetic braking and gravitational radiation. 

Key words, stars: neutron - stars: individual (X1822-371) — X-rays: binaries — X-rays: pulsars 
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£S i 1- Introduction in 2000 (Obs ID: 671) and in 2008 (Obs ID: 9076 and 9858), 

, already included in the work of BU10, estimated a value of the 

Cn . X1822-371 is an eclipsing compact binary system with a period orbkal period derivative of (0 83 ± 16) x 10 -io s/s> with the 

■ of 5.57 hr hosting a 0.59 s X-ray pulsar. Several authors have errQr at the gQ% confidence level almost a factor of two smaller 

j reported new orbital ephemeris of th e source using observ ations man ^ varae re p 0r t e d by BU10 

(N I performed in different energy bands. Bur deri et al. | feM here- We summarise the values of the ecli re fe re nce time T: 

after BU 10) analysed X-ray data of XI 822-371 covering the pe- , ,. , . , „ , , ,. , . , , . . • 

■ • riod from ioo 6 t0 ?008 to determine the eclicse times of the the orbltal penod P ° lb0 ' and the orbltal penod denvatlve P < 
^ , nou rrom iyyo to zuuo io determine me ecnpse times or tne . mnn n a in „„a m 1 ;„ XoKi^m 



orb 

, . , x , obtained by PA00, BUI 0, BA 10, and Jill in TableQ] 

• »-h source and improved the previous X-ray ephemeris of XI 822- , ,. , , • , • , , ^ 

><] " 371 reported by Parmar et al. (2000, hereafter PA00) that cov- } n thls , ,, w ° rk ' , we comprehensively examine both X-ray 

U i ered the period from 1977 to 1996. BU10 added their data to and ~kp« "tnyal ^ UP 

*u j u mnAt j- •<■• a ■ t - t.i, u-* i ephemeris of XI 822-371, adding to the eclipse arrival times re- 

, . w t those used by PA00 finding a positive derivative of the orbital r ., „,,,„. . f. _ \ 

period of (1.499 ± 0.071) x 10" 10 s/s that is compatible with the P™**"? B ™ ™ ° n f ?^ ? * S , observatlon P e i 1 "- 

previous one given by PA00 but with a smaller associated error. formed ln 20 f 06 - W ^ . ds ° ™ lud f a data P 0lnt f mm a ob " 

17^ — — rn j ^Xml i. nnm , , servation performed in 1989, and a data point from a ROS AT ob- 

Bavlessetal. (2010, hereafter BA10) obtained the opti- e . ,. , r . nr , , ,. 

iJTTTt — T~ ■ Tvion t-71 i . , servation performed in 1992. We critically examine the discrep- 

cal/UV ephemeris of XI 822-371 using data covering the period • , , ,. , , • ■ , ■ , , • • 

c mnn < onn/c • , , c ?•. , . , ancies that have emerged in calculating the orbital ephemeris in 

from 1979 to 2006. They obtained a value of the orbital period . , ? „ , , , • ,- , 



previous papers and, finally, show the ephemeris of the X1822- 
371 by combining the optical/UV and X-ray data-sets. 



derivative of (2.12 ± 0.18) x 10 10 s/s, which is compatible with 
that reported by PA00 but slightly larger than the value proposed 
by BUlO. 

Ui et alJ (1201 lL hereafter JI 1 1 ), using the X-ray eclipse arrival 
times reported by PA00 and the eclipse arrival times inferred by 
the two Chandra/HETG observations of XI 822-371 performed Suzaku observed XI 822-371 on 2006 October 2 with an 
elapsed time of 88 ks, the start and stop times of the obser- 



2. Suzaku observation 



Send offprint requests to: R. Iaria, e-mail: iaria@fisica.unipa. it vation corresponding to 54010.48 and 5401 1.50 MJD, respec- 
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Table 1. Journal of the ephemerides of X1822-371 discussed in this work. 



Parameters 


Parmar et al. (2000) 


Burderi et al. (2010) 


Bayless et al. (2010) 


Jietal. (2011) 


n (mjd q ) 


45614.80964(15) 


45614.80948(14) 


45614.81166(74) 


45614.80927(25) 


PorhQ (S) 


20054.1990(43) 


20054.2056(22) 


20054.1866(69) 


20054.2181(41) 


PoaCxIO-'Os/s) 


1.78(20) 


1.499(71) 


2.12(18) 


0.827(95) 


X 2 /(d.o.f.) 


21.4/16 


38.69/25 


70.04/32 


35.99/19 



Note — Uncertainties are at the 68% c. 1. for a single parameter. We show the reference time Ti of the eclipse arrival times in units of MJD, the 
orbital period f or bo in units of seconds calculated at Tg, the derivative of the orbital period P or \, in units of s/s, and finally the % 2 /(d.o.f.) obtained 
fitting the eclipse arrival times with a quadratic function. 



tively. Both the X-r ay Imaging Spectrometer (0.2-12 keV, XIS; 
iKovama et ai]|2007l) and the Hard X-ray Detector (10-600 keV, 
HXD: lTakahashi et al.l2007l) instruments were used during these 
observations. In this work, we used only the XIS data. There are 
four XIS detectors, numbered to 3. The XIS0, XIS2, and XIS3 
detectors use front-illuminated CCDs and have very similar re- 
sponses, while XIS1 uses a back-illuminated CCD. 

We reprocessed the observation using the aepipeline tool in- 
cluded in the Suzaku FTOOLS Version 16 applying the latest 
calibration available as of 2011 March. During the observation, 
XIS0 and XIS 1 were used adopting the quarter window option 
(frame time 2 s), while XIS2 and XIS3 worked in full window 
(frame time 8 s) mode. We barycentred the XIS data using the 
SUZAKU tool aebarycen and adopting as the best estimate of the 
source coordinates those derived from the 2008 Chandra obser- 
vations (RA: 18 25 46.81, DEC: -37 06 18.5, uncertainty: 0.6"). 

We extracted the four XIS light curves in the 1-10 keV en- 
ergy band selecting a circular region centred on the source. We 
adopted a radius of 130 pixels for XIS0 and XIS1 and 160 pix- 
els for XIS2 and XIS3. The four light curves are quite similar 
and enclose four orbital periods of XI 822-371, thus we used the 
FTOOL lcmath to combine the four XIS light curves. The com- 
bined XIS light curve is shown in Fig. Q] adopting a bin time of 
128 s. 



Bin time: 128.0 s 
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Fig. 1. Combined XIS light curve of X1822-371 in the 1-10 keV 
energy band. The adopted bin time is 128 s. 



Table 2. Journal of all the available X-ray eclipse times. 



P /"* 1 1 r\ c & ' 1 i tn & ( T Pi i 
J_/dipaC 11111C \JLJq) 


Error 
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Q/J 8 A 
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tlmU- 1 oCdll 
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Einstein 


1AAA\ ^1 


U.UUjU 






Einstein 


2445580.4932 


0.0005 


-150 


j 


EXOSAT 


2445615.30940 


0.00038 







EXOSAT 


2445963.00914 


0.00033 


1498 




EXOSAT 


2445963.24046 


0.00030 


1499 




EXOSAT 


2445963.47254 


0.00034 


1500 




EXOSAT 


2446191.63643 


0.00031 


2483 




EXOSAT 


2446191.86768 


0.00033 


2484 




EXOSAT 


2446192.10008 


0.00029 


2485 




EXOSAT 


2447760.22900 


0.00030 


9241 




Ginga 


2448692.84396 


0.00070 


13259 




ROSAT 


2449268.00984 


0.00040 


15737 


3 


ASCA 


2450353.35425 


0.00035 


20413 


3 


ASCA 


2450353.58728 


0.00023 


20414 


3 


RXTE 


2450701.51870 


0.00120 


21913 


3 


BeppoSAX 


2450992.58580 


0.00230 


23167 


4 


RXTE 


2451780.13170 


0.00190 


26560 


4 


Chandra 


2451975.56934 


0.00056 


27402 


4 


XMM-Newton 


2451975.56935 


0.00031 


27402 


4 


RXTE 


2452432.59458 


0.00030 


29371 


4 


RXTE 


2452488.53300 


0.00038 


29612 


4 


RXTE 


2452519.63569 


0.00085 


29746 


4 


RXTE 


2452882.65470 


0.00037 


31310 


4 


RXTE 


2454011.17300 


0.00090 


36172 


5 


Suzaku 


2454607.69592 


0.00056 


38742 


4 


Chandra 



Note — References. (1) iHellier & Smaldll994 (2 ) iHellier & Masonl 
119891, (3) lParmar et al.ll200Cl, (4) lBurderi et alJl2010L (5) this work. The 
number of cycles for each eclipse time is discussed in sec. 13. II 



During only the third orbital passage of XI 822-371, the 
eclipse was fully covered by Suzaku at a time of 55,000 s from 
the start time. To estimate the eclipse arrival time, we folded 
the combined XIS light curve, adopting the ephemeris reported 
by BU10 and a bin time of 128 s. We fitted the orbital light 
curve to derive eclipse arrival times by adopting the same pro- 
cedure described in BU10, obtaining an eclipse time passage at 
54, 010.6730 + 0.0009 MJD with an associated error at the 68% 
confidence level. 

3. The ephemeris of X1 822-371 

For clarity's sake, we show in Table [2] the X-ray eclipse arrival 
times that we used to update the X-ray ephemeris of XI 822-371. 
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Most of these data points were included in the timing analysis 
of BU10. To their data set, we added eclipse arrival times from 
Ginga (1989), ROSAT (1992), and, most importantly, Suzaku 
(2006). The ephemerides showed in Table Q] and in the anal- 
ysis now described are given in barycentric dynamical time. 
We note that the RXTE arrival times from 1998 to 2003 re- 
ported in Tab. 1 of BU10 (except for the second point corre- 
sponding to cycle 23,167) are not the eclipse arrival times, as 
erroneously stated, but the times of passage through the as- 
cending node (which differs from the eclipse time by P rb/4). 
Nevertheless, the corresponding RXTE time delays were cor- 
rectly shown in Fig. 1 of BU10 and correctly used to derive the 
orbital ephemeris, which are therefore unaffected by this mis- 
take. The correct RXTE eclipse arrival times are shown in our 
Table |2] The X-ray ephemeris of XI 822-371 reported by BUI 
and Jill show a large discrepancy in the quadratic term by al- 
most a factor of two. JI1 1 suggested that the discrepancy in the 
time delay associated with the last two Chandra observations is 
caused by BU10 not folding the Chandra light curves to estimate 
the eclipse arrival time, which instead was done by BU10. To un- 
derstand the reason for this discrepancy, as a first step we tried to 
reproduce the results of JI1 1 by using the same data they used in 
their analysis. These consist of a total of 22 eclips e -times , which 
are, respectively, those given bv lHellier & Smalel ( fl994b . PA00, 
and three eclipse arrival times obtained from three Chandra ob- 
servations corresponding to obsID 671, 9076, and 9058 derived 
by JI1 1 (see Tab. 2 in their paper). 

We found the corresponding time delays following their pro- 
cedure, namely we determined the ti me delays with re s pect t o 
the best-fit linear ephemeris shown by iHellier & Smald d 19941) . 
that is 

= 2445615.30942(14)JD Q + 0.232109017(33)^, 

and fitted the time delays with a quadratic function. We ob- 
tained best-fit values consistent with the ones reported in Jill. 
We showed in Fig. [2] the time delays in unit s of days associ- 
ated w ith the eclipse arrival times shown by IHellier & Smald 
i 19941) and PA00 with black open squares. The delay times as- 
sociated with the Chandra eclipse arrival times showed by JI1 1 
were plotted using red diamonds. The dashed line corresponds 
to the quadratic best-fit curve given by JI1 1 . 

After establishing the reproducibility of the parameter's es- 
timates of Jill, we explored the cause of the discrepancy in 
the fitting results given by BU10 extracting the eclipse arrival 
times from each Chandra barycentred and folded light curve (ob- 
sIDs 671, 9076, and 9858). We show the Chandra eclipse arrival 
times in Table [3] The corresponding delays were estimated as 
described above. 

Table 3. Eclipse arrival times of the three Chandra observations. 



ObsID 


Eclipse Time (MJD Q ) 


Cycle 


Delay (days) 


671 


51779.8638(11) 


26561 


0.0067(11) 


9076 


54607.19610(56) 


38742 


0.01914(56) 


9858 


54609.74890(31) 


38753 


0.01875(31) 



Note — Column 1: Chandra obsIDs. Column 2: our estimation of the 
eclipse arrival time. Column 3 : the corresponding cycl e with respect to 
the linear ephemeris given bv lHellier & Smald (1994). Column 4: the 
corresponding delay in units of days. The errors are at the 68% c.l. 

We found that the discrepancies of the eclipse arrival times 
between our analysis and JI1 l's are -40 ± 140 s, 440 + 70 s, and 



in rr 




-10 4 1 4 2xl0 4 3x10" 

N (cycles) 

Fig. 2. Time delays of th e eclipse arrival t imes w ith respect 
to the linear ephemeris of IHellier & Smalel dl994l) . The black 
open squares cor r espon d to the eclipse arrival times shown by 
IHellier & Smald (11994 and PA00, the red diamonds to the 
Chandra eclipse arrival times given by Jill, the green filled 
squares to the Chandra eclipse arrival times reported in this pa- 
per (see Table [3), the blue open circles to the Chandra eclipse 
arrival times reported by BU10. The dashed and solid lines are 
the best-fit quadratic curve obtained by Jill and in this paper 
(see text), respectively. 

440+50 s for obsIDs 671, 9076, and 9858, respectively. Since the 
errors are at the 68% c. 1., the eclipse arrival times correspond- 
ing to the obsIDs 9076 and 9858 are not compatible. We show 
our Chandra delays with green filled squares in Fig. [2] Only two 
eclipse times were derived by BU10 from the Chandra obser- 
vations, because light curves of obsID 9076 and obsID 9058 
were combined to obtain a single folded light-curve and a sin- 
gle eclipse-time passage, with a smaller uncertainty, since the 
observations were sufficiently close in time to each other. We 
show the two corresponding time delays with blue open circles 
in Fig. [2] We note that our delays and those given by BU10 are 
widely compatible. 

Fitting the t ime delays correspond ing to the eclipse arrival 
times given by IHellier & Smald d 19941) . PA00, and our three 
eclipse arrival times reported in Table [3] for a total of 22 data 
points, we obtained 

T ecl = 45614.80954(14)MJD o + 0.2321088628(21)JV 

+ 1.648(72) x lO^'JV 2 , (1) 

with a x 2 /(d.o.f.) = 25.6/19 and the errors are at the 68% 
confidence level, the uncertainties in the parameters have been 

scaled by a factor ^jx^ ed to take into account a x^ ed of the best- 
fit model larger than 1 . We note that the quadratic term is larger 
than that shown by Jill. The corresponding P O rb0 and P OI \, are 
20054.20574(17) s and 1. 420(63) x 10~ 10 s/s, respectively. These 
values are compatible within one <x with the ones given by BU10 
(see Table [TJ. 

3. 1 . Updated X-ray ephemeris of X1 822-371 

As a first step, we found the X-ray ephemeris of XI 822-371 
using the eclipse arrival times adopted by BU10 excluding the 
Chandra eclipse arrival times and including the eclipse arrival 
times taken with Ginga, ROSAT, and Suzaku (see Tab. 13 for a 
total of 28 available data points. The Suzaku data-point is very 
important in this respect, since it was taken in 2006 and there- 
fore precedes the last Chandra data-points taken in 2008. This 
is very important to fill the time gap between the last RXTE ar- 
rival time taken in 2003 and the most recent Chandra observation 
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taken in 2008, and therefore gives us the opportunity to discrimi- 
nate more clearly between the Chandra eclipse arrival time as re- 
ported by JI1 1 and our measurement (which is compatible with 
the one reported by BU10). 

We found the delays of the eclipse arrival times by subtract- 
ing from our measurements the eclipse arrival times predicted 
by a constant orbital period model adopting the orbital period, 
P or b o, and the reference time, Tf., given by PA00. The time de- 
lays were plotted versus the orbital cycle number N. The integer 
N is the exact number of orbital cycles elapsed since TS; the cy- 
cle number N corresponding to each eclipse arrival time is shown 
in column 3 of Table [2] We then fitted the time delays using a 
parabolic function obtaining ax 2 /(d.o.f.) of 33.63/25. We found 
that T e Q = 45614.80959(16) MJD , P olb0 = 20054.2020(28) s, 
and P or b = 1.626(90) x 10' 10 s/s with the associated errors at the 
68% confidence level. All these values are compatible within 
one cr with those given by BU10 (see column 3 in Table [TJ, this 
suggests that the Chandra eclipse arrival times given by BU10 
are in agreement with all the previous points. 

To update the X-ray ephemeris of XI 822-371, we then in- 
cluded the Chandra eclipse arrival times given by BU10 in our 
data set for a total of 30 available data points. We found the cor- 
responding delays and cycle numbers as described above. The 




' 2xl0 4 

N (cycles) 

Fig. 3. Upper panel: Eclipse time delays with respect to a con- 
stant orbital period model plotted versus the orbital cycle for all 
the available X-ray eclipse time measures together with the best- 
fit parabola. Lower panel: residuals in units of cr with respect 
to the best-fit parabola. The black full squares points are from 
BU10, the red diamonds are the data added in this work. 

time delays are shown in the upper panel of Fig. [3] We plotted 
the time delays used by BU10 with black full squares, while the 
time delays added in this work and corresponding to the Ginga, 
ROSAT, and Suzaku eclipse times are shown with red diamonds. 
We then fitted the time delays using a parabolic function, result- 
ing in the ephemeris 

T ed = 45614.80953(16)MJD o + 0.232108853(30)/V 

+ 1.757(93) x 10" 1 l N 2 , (2) 

where the associated errors are at the 68% confidence level. We 
obtained ay 2 1 (d.o ./.) of 4 1 .2/27, the best-fit curve is shown with 
a solid line in Fig. [3] We show the residuals in units of cr in the 
lower panel of Fig. [3] and report the obtained values of TZ, P 01 -b o, 
and P or b in the second column of Table [4] 



We found that the derivative of the orbital period, P or b, is 
1.514(80) x 10 -11 s/s, compatible with the value of 1.499(71) x 
10~' 1 s/s estimated by BU10. 



Table 4. Updated X-ray and optical/UV ephemeris of XI 822- 
371. 



Parameters 


X-ray 


Optical/UV 


T e (MJD ) 


45614.80953(16) 


45614.8116(11) 


^orb (S) 


20054.2049(26) 


20054.188(10) 


Arb 


1.514(80) 


2.10(26) 


X 2 1 (d.o.f.) 


41.2/27 


71.38/33 



Note — Uncertainties are at the 68% c. 1. for a single parameter. The 
parameters are defined as in Table [T] the derivative of the orbital period 
is in units of 10~ 10 s/s. The updated values of Tg, P OT bo, an d Ant using 
X-ray data and optical/UV data are shown in Cols. 2 and 3, respectively. 



3.2. Updated optical/UV ephemeris of X1 822-371 

BA10 used 35 optical/UV eclipse arrival times shown in Tab. 1 
of their paper to find the best-fit optical/UV ephemeris of X 1 822- 
37 1 given by 

T ed = 45614.81 166(74)MJD + 0.2321 0864 1(80)/V 

+ 2.46(21) x 10~ n /V 2 , (3) 

where the errors are at the 68% confidence level (Bay less, private 
communication). We added to their data the optical eclipse ar- 
rival ti me 2,443,629.841 +0.013 JD Q given bv lHellier & Masonl 
( 1 989) and not included in BA 1 0. 

Using the 36 optical/UV data points and following the proce- 
dure described in the previous section we found the correspond- 
ing time delays. Fitting them with a parabola, we obtained the 
following optical/UV ephemeris 

T ecl = 45614.81 16(1 1)MJD + 0.23210865(12)/V 

+ 2.44(31) x 10~ n /V 2 , (4) 

with a x 2 1 (d.o.f.) of 71.38/33 and the errors are at the 68% 
confidence level. The uncertainties have been scaled by a fac- 
tor Jx^ed to take mt0 account a x 2 ec i °f tne best-fit model larger 
than 1. This explains why the uncertainties in the optical/UV 
ephemeris we have shown are larger than the ephemeris shown 
by BA10. The updated optical/UV ephemeris are consistent with 
those given by BA10. We report the corresponding values of T e Q , 
P 0I b o, and P m \, in the third column of Table|4] In the upper panel 
of Fig. |U we show the time delays for each eclipse arrival time 
of XI 822-371 for the X-ray (red full squares) and optical/UV 
bands (black full squares) for a total of 66 data points. The solid 
and dashed lines correspond to the best-fit parabolas reproduc- 
ing the X-ray and optical/UV ephemerides showed in eqs.|2]and 
|U respectively. 

We compare the X-ray and optical/UV residuals with the 
X-ray best-fit parabola in Fig. [4] (middle panel). Although we 
plot the residuals for the best-fit parabola obtained from the 
X-ray time delays, we note that almost all of the optical/UV 
data are close to the best-fit curve. The largest discrepancies 
are associated wi t h the last two optical eclipse times shown by 
iHellier & Masonl (1 19891) corresponding to orbital cycles 7,243 
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Fig. 4. Top panel: the optical/UV (black filled squares) and 
X-ray (red filled squares) time delays. The dashed and solid 
lines correspond to the optical/UV and X-ray best-fit parabolic 
curve. Middle panel: residuals with respect to the X-ray best- 
fit parabolic curve. Bottom panel: residuals with respect to the 
optical/UV best-fit parabolic curve. 



and 7,600 and the two UV eclipse arrival times obtained with 
HST and reported by BA10; these last two data points are at 
orbital cycles 35,387 and 35,395, respectively. All the other op- 
tical/UV points are within two cr of the corresponding values of 
the best-fit X-ray ephemeris. 

In Fig. 3] (bottom panel), we show the X-ray and optical/UV 
residuals with respect to the optical/UV best-fit curve. In this 
case, the X-ray data are mainly below the best-fit optical/UV 
parabola. 

3.3. Time-lag between optical/UV and X-ray eclipse times 

To our knowledge, there are only two simultaneous X-ray and 
optical observations of the eclipse of XI 822-3 71 reported by 
IHellier & Masonl d!989h and IHellier etafl d!990l) : these authors 
showed that the optical eclipse times lag the X-ray eclipse times 
by 3.0 + 3.4 min, and 180 + 50 s, respectively. 

The optical eclipses are also wider than the X-ray eclipses; 
the different width suggests a different origin f o r the o ptical and 
X-ray eclipses, respectively. Hellier & Mason (1989) proposed 
that the X-ray emission comes from an accretion disc corona 
(ADC) with a radius half of the outer accretion disc radius, while 
the optical emission is produced by a more extended disk struc- 
ture. Furthermore, the optical eclipse lags the X-ray eclipse be- 
cause of the asymmetric disk structure probabl y caused by the 
stream impact onto the outer accretion disk. IHellier & Masonl 
(fl989l) . modelling the X-ray and optical light curves of XI 822- 
371, found an optical eclipse time-lag of ~ 0.01 in units of 
orbital phase, corresponding to a time-lag of 200 s. BA10 dis- 
cussed a marginally significant time-lag between the optical/UV 
and X-ray ephemeris of 100 + 65 s and 122 s with respect to the 
X-ray ephemeris reported by PA00 and BU10, respectively. 

Since we have used an unprecedentedly large amount of op- 
tical/UV and X-ray eclipse times, we can estimate the average 



time-lag along 50,000 orbital cycles with good accuracy. We 
fitted simultaneously the X-ray and optical/UV time delays al- 
lowing the constant terms of each parabola free to vary and 
constraining the values of the linear and quadratic parameters 
of each parabola to the same value, since the orbital period of 
XI 822-371 and its derivative cannot depend on the considered 
waveband. 

Fitting the time delays, we obtained a large x 2 /(d.o.f.) of 
124.04/62 and found that the best-fit values of the linear and 
quadratic terms are (3.7 + 2.8) x 10~ 3 s and (1.595 + 0.086) x 10~ 6 
s, respectively. The constant terms are 121 + 36 s and -6 + 16 s 
for the optical/UV and X-ray data-sets, respectively. Using these 
values, we obtained the ephemerides for the X-ray and opti- 
cal/UV data 



+ 1.847(99) x 1Q- U N 2 , 



T edeptm = 45614.81 104(42)MJD o + 0.232108828(32)JV 



(5) 



+ 1.847(99) x 10" 1 l N 2 . (6) 

The corresponding orbital period derivative is 1.591(86) x 
10~ 10 s/s, and the reference time T is 45614.80957(19) and 
45614.81 104(42) MJD for X-ray and optical/UV data-sets, re- 
spectively; all the errors are at the 68% c.l. For clarity's sake, 
the values of Tf., P i-bo, and P or b are showed in Table [5] In the 
upper panel of Fig. [5] we show the X-ray (red points) and the op- 
tical/UV (black points) time delays; the dashed and solid curves 
are the optical/UV and X-ray best-fit parabolas, respectively. 

From our analysis, we found a time-lag of 127 + 52 s, which 
is significant at a confidence level of 2.4 cr. In the bottom panel 
of Fig. [5] we show the residuals of the X-ray (red points) and 
optical/UV (black points) delays with respect to the X-ray best- 
fit parabola. 




1 tiitI itii IittiItittItiii 


1 1 1 1 1 1 1 1 1 1 1 1 1 

'^fj til ljt t ; 


1 1 1 1 1 1 1 1 1 1 1 1 

1 ■ 


-1 , ^ , , , 1 , , 





-1() 4 



lu 4 2xl0 4 
Orbital cycles 



3xl0 4 



Fig. 5. Upper panel: the optical/UV (black filled squares) and 
X-ray (red filled squares) time delays fitted with two parabo- 
las having the same linear and quadratic terms. The solid and 
dashed parabolas correspond to the X-ray and optical/UV best- 
fit curves. Lower panel: Residuals in units of cr with respect to 
the best-fit parabola describing the X-ray ephemeris shown in 
eq.0 



The values of the optical/UV time-lags with respect to the 
best-fit parabola describing the X-ray ephemeris given by eq. [5] 



6 R. Iaria et al.: Detailed study of the orbital ephemeris of XI 822-371 

Table 5. Ephemeris of XI 822-371 fitting simultaneously X-ray and optical/UV data. 



1 Ox-ray 


je 

%pticalfUV 


-^orbO 




X 2 /(d.o.f.) 


(MJD ) 


(MJD Q ) 


(s) 


(xlO^ 10 ) s/s 




45614.80957(19) 


45614.81104(42) 


20054.2027(28) 


1.591(86) 


124.04/62 



Note — Uncertainties are at 68% c. 1. for a single parameter. The parameters are defined as in Table [Tj The uncertainties in the parameters have 
been scaled by a factor Jx^ ed to take into account &X^ ed of the best-fit model larger than 1. 



are shown in Table [6] We note that the largest optical time-lags 
are associated wi t h the last two optical eclipse times shown by 
iHellier & Mason! d!989l) corresponding to orbital cycles 7,243 
and 7,600 and with the recent two UV eclipse arrival times ob- 
tained with HST and reported by BA10; the corresponding time- 
lags are 433 + 86 s, 282 + 86 s, 364 + 86 s, and 452 + 86 s, 
respectively. These va lues are larger than 200 s, in disagreement 
with that predicted bv lHellier & Mason] d!989l) modelling the X- 
ray and optical light curves of XI 822-371. 




2000 4000 6000 8000 

t-2,443,629. 8410 (days) 



Fig. 6. Upper panel: the optical/UV (red filled squares) time- 
lags fitted with a constant (dashed line) and a sinusoidal func- 
tion f(t) (solid line) with a period of 18 yrs (see text). Lower 
panel: Residuals in units of <x with respect to the best-fit sinu- 
soidal function. 



Table 7. Best-fit parameters of the sinusoidal modulation fitting 
the optical/UV time-lags. 



Parameters 


best-fit 1 


best-fit 2 


A(s) 


161 ±24 


105 ± 24 


B(s) 


194 ± 25 


-267 ± 43 


P (days) 


6593 ± 452 


283.1 ±0.6 


t (days) 


-1180 ± 481 


239 ± 16 


X 2 (Aof) 


34.52(32) 


35.57(32) 



Note — Uncertainties are at the 68% c. 1. for a single parameter. In 
the second and third columns, we show the best-fit parameters for the 
modulation of 18 years and 283 days, respectively. 



We fitted the optical/UV time-lags with a constant obtain- 
ing a large x 2 (d.o.f.) of 81.25(35); the constant value was 127 s, 
which is similar to the averaged time-lag previously discussed. 
In Fig. [6] (top panel), we show the optical/UV time-lags as a 
function of time in units of days, the dashed line being the con- 
stant function. We note that if we remove the four optical/UV 
points mentioned above, we find no significant time-lag between 
the optical/UV and X-ray eclipse times, with a best-fit value of 
31+ 46 s. 

Because we found a large value ofx 2 red , after a visual inspec- 
tion of the fit residuals, we decided to fit the time-lags with the 
function fit) — A - B sin [2n/P(t - to)]- In this case, we largely 
improved the fit for two different sets of parameters. For the 
first set, we obtained ax 2 (d.o.f.) = 34.52(32) and a probabil- 
ity of chance improvement with respect to the fit with a con- 
stant of 4.03 x 10~ 6 . The values of the best-fit parameters are 
A = 161+24s,B = 194±29 s,P = 6593 ± 452 d (18.1 ± 1.2 yr), 
and t = -1 180 + 481 d, the errors being at the 68% c.l.. For the 
second set of parameters, we obtained a x 2 (d.o.f.) = 35.57(32) 
and a probability of chance improvement with respect to the fit 
with a constant of 6.45 x 10~ 6 ; in this case, A = 105 ± 24 s, 
B = -267 ± 43 s, P = 283.1 + 0.6 d, and t = 239 + 16 d, the 
errors being at the 68% c.l.. The best-fit values of both the fits 
are shown in Table Q 



o 




2000 4000 6000 



t-2,443,629.8410 (days) 

Fig. 7. Upper panel: the optical/UV (red filled squares) time-lags 
fitted with a sinusoidal function /(f) (solid line) with a period of 
283 days (see text). Lower panel: Residuals in units of <x with 
respect to the best-fit sinusoidal function. 
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Table 6. Journal of the optical/UV time-lags 
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Eclipse Time 


Time-lag 


residuals best-fit 1 


residuals best-fit 2 


(JD Q -2,400,000) 


(s) 


(s) 


(s) 


43629.8410 + 0.0130 


-951 ± 1123 


-938 ± 1123 


-1278 ± 1123 


44044.8450 ± 0.0060 


-1477 ±518 


-1444 ±518 


-1395 ±518 


44090.1140 + 0.0080 


-800 ± 691 


-767 ± 691 


-644 ± 691 


44101.0280 ±0.0080 


-377 ±691 


-344 ± 691 


-241 ± 691 


44105.6650 ±0.0060 


-824 ±518 


-791 ±518 


-701 ±518 


44106.5970 ± 0.0060 


-516 ± 518 


-483 ± 518 


-396 ±518 


44137.9350 ±0.0100 


-227 ± 864 


-195 ±864 


-252 ± 864 


44411.1320 ±0.0080 


218 ±691 


242 ± 691 


247 ± 691 


44412.0580 ±0.0080 


8 ±691 


32 ± 691 


32 ± 691 


44664.8120 ±0.0080 


-1057 ±691 


-1051 ±691 


-916 ±691 


44783.8940 ± 0.0060 


-172 ± 518 


-179 ±518 


-542 ±518 


45579.5650 ± 0.0030 


16 ±259 


-114 ±259 


-163 ±259 


45580.7250 ± 0.0030 


-31 ±259 


-162 ±259 


-216 ±259 


45615.3115 ±0.0020 


166 ± 173 


-30 ± 173 


-173 ± 173 


45937.7110 ±0.0021 


192 ± 181 


-4± 181 


-160 ± 181 


46234.5787 ±0.0018 


228 ± 156 


-20 ± 156 


-83 ± 156 


47296.4798 ± 0.0010 


433 ± 86 


80 ±86 


158 ± 86 


47379.3410 ± 0.0010 


282 ± 86 


-72 ± 86 


25 ±86 


47999.7674 ± 0.0039 


163 ± 337 


-160 ±337 


209 ± 337 


49163.0960 ±0.0028 


-114 ±242 


-253 ± 242 


37 ± 242 


49164.0237 ±0.0028 


-178 ±242 


-317 ±242 


-25 ± 242 


49164.2542 ±0.0028 


-317 ±242 


-456 ± 242 


-164 ±242 


49165.1852 ±0.0028 


-96 ± 242 


-234 ± 242 


59 ± 242 


49166.1190 ±0.0028 


368 ± 242 


229 ± 242 


524 ± 242 


50250.5308 ± 0.0008 


49 ±69 


65 ±69 


22 ±69 


50250.7626 ± 0.0008 


22 ±69 


38 ±69 


-3 ±69 


50252.6196 ±0.0008 


33 ±69 


49 ±69 


18 ±69 


51264.8446 ± 0.0035 


-413 ± 302 


-409 ± 302 


-713 ± 302 


51373.7094 ±0.0008 


47 ±69 


40 ±69 


-32 ± 69 


51756.4577 ±0.0012 


-15 ± 104 


-74 ± 104 


78 ± 104 


51782.4542 ±0.0012 


3± 104 


-60 ± 104 


-34 ± 104 


52089.7692 ± 0.0008 


139 ± 69 


24 ±69 


-40 ± 69 


53618.6767 ±0.0050 


58 ± 432 


-285 ± 432 


-134 ±432 


53828.9720 ± 0.0010 


364 ± 86 


11 ±86 


11 ±86 


53830.8299 ±0.0010 


452 ± 86 


99 ±86 


95 ±86 


53932.7201 ±0.0040 


-80 ± 346 


-434 ± 346 


-93 ± 346 



Note — The optical/UV eclipse times (1st column), the corresponding time-lags (2nd column) with respect to the best-fit parabola describing the 
X-ray ephemeris showed in eq.[5] the time-lags after removal of the sinusoidal modulations (3rd and 4th column, respectively). The errors are at 
the68%c.l. 



The best-fit curve corresponding to the first set of parame- 
ters (hereafter best-fit 1) is shown in Fig. [6] (top panel) with a 
solid line. In the bottom panel of Fig. [6j we show the residu- 
als with respect to the sinusoidal modulation in units of cr. The 
values of the residuals corresponding to the first set of parame- 
ters are shown in the third column of Table [6] The best-fit curve 
corresponding to the second set of parameters (hereafter best-fit 
2) is shown in Fig. [7] (top panel) with a solid line. In the bot- 
tom panel of Fig. |7J we show the residuals with respect to the 
sinusoidal modulation in units of cr. The values of the residu- 
als corresponding to the second set of parameters are shown in 
the fourth column of Table [6] We checked our results fitting the 
X-ray time-lags with respect to the best-fit parabola giving the 
X-ray ephemeris. In this case, the residuals do not show any si- 
nusoidal modulation, as expected. Fitting them with a constant, 
we found & x 2 {d.o.f.) = 42.8(29). 

The two best-fit sinusoidal functions indicate that the opti- 
cal/UV eclipse lags the X-ray eclipse with a time shift of either 
161 ± 24 s or 105 + 24, which are significant at a confidence level 
of 6.1 cr and 4.4cr, respectively. These values are compatible with 



the time-lag predicted bv lHellier & Masonl d!989l) . of ~ 200 s. 
Our most intriguing result, never previously detected, is that the 
optical/UV eclipse times may oscillate in time with an amplitude 
of either 194 + 29 s (best-fit 1) or 267 + 43 s (best-fit 2), these 
values being significant at a confidence level of 6.7<x and 6.2cr, 
respectively. The detected periods are ~ 18 yr (best-fit 1) and 
~ 283 days (best-fit 2). Both the detected periods are very long 
and d ifficult to explain by invok ing a superhump phenomenon 
(see I Wang & ChakrabartwOlOl and reference therein). The su- 
perhump excess e is defined as P s h/P rb - 1, whe re P s h is the su- 
perhu mp period, and also as e = 0.1 8g + 0.29g 2 (Patterson et al. 
2005), where q is the mass ratio milm\ with m\ the neutron-star 
mass. Since the mas s function and inclination an gle of X1822- 
371 is well known dJonker & van der Klis 2001), assuming a 
neutron-star mass of 1 .4 M Q and an inclination angle of 87° we 
find that q =* 0.29 and e ~ 0.077. Consequently the superhump 
period should be = \.QllP or b = 21598.38 s. A possible beat 
phenomenon between the superhump period and the orbital pe- 
riod could produce a period given by (l/P or b - l/Psh)' 1 - 3.25 
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days. This value is shorter than the two periodicities that we re- 
ported. 

We found that the lag changes with time and that this varia- 
tion is compatible with a sinusoidal modulation at two different 
periods of 18 yr and 283 days, respectively. However, we cannot 
exclude shorter periodicities, but our data set of only 36 opti- 
cal/UV eclipse times spanning a time of 10,000 days do not al- 
low a rigorous study. To investigate this aspect, long optical/UV 
observations of XI 822-371 covering several contiguous orbital 
periods of the source would be necessary. 

Finally we note that both these best-fit curves are strongly 
driven by the optical measures corresponding to orbital cycles 
7,243 and 7,600 and by the recent two UV eclipse arrival times 
obtained with HST and reported by BA10; these last two points 
are at orbital cycle 35,387 and 35,395. 
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4. Conclusions 

We have revisited and discussed the X-ray and optical/UV 
ephemerides of XI 822-371. Fitting simultaneously the opti- 
cal/UV and X-ray time delays, we have found that the opti- 
cal/UV eclipses of X 1 822-37 1 lag the X-ray eclipses by 1 27 + 52 
s with a significance level of 2.4 <x. However, this time-lag may 
not be constant in time. Fitting the optical/UV time-lags, we 
have found a statistically significant variation, which is compat- 
ible with a sinusoidal modulation at two different periods, ~ 18 
yr and 239 d. In the first case, the optical/UV eclipses lag the X- 
ray eclipse by an average time of 161 s (significance 6.7cr) and 
this delay oscillates in time around this value with an amplitude 
of 194 s (significance 6.7<x). In the second case, the optical/UV 
eclipses lag the X-ray eclipse by an average time of 105 s (signif- 
icance 4 .4cr), and this delay oscillates in time around this value 
with an amplitude of 267 s (significance 6.2cr). 

Owing to the relatively small number of points over a long- 
time span of 30 yr, we cannot be sure of the period of this mod- 
ulation, because we cannot exclude much shorter periods. Long 
and relatively continuous optical/UV observations are necessary 
to prove or disprove the presence of this periodicity in the opti- 
cal/UV eclipse time-lags. 

Our results confirm the value of the orbital solution derived 
by the X-ray eclipse times given by BU 1 and that the orbital pe- 
riod derivative is three orders of magnitude larger than expected 
on the basis of the conservative mass transfer driven by magnetic 
braking and gravitational radiation. We have also confirmed this 
result by combining the X-ray data and the optical/UV data of 
XI 822-371. 
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